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We study a weighted-set graph coloring problem in which one assigns q colors to 
the vertices of a graph such that adjacent vertices have different colors, with a vertex 
weighting w that either disfavors or favors a given subset of s colors contained in 
the set of q colors. We construct and analyze a weighted-set chromatic polynomial 
Ph(G,q,s,w) associated with this coloring. General properties of this weighted- 
set chromatic polynomial are proved, and illustrative calculations are presented for 
various families of graphs. This study extends a previous one for the case s = 1 and 
reveals a number of interesting new features. 



I. INTRODUCTION 

Recently, two weighted graph coloring problems have been formulated and studied in 
which one assigns q colors to the vertices of a graph subject to the condition that adjacent 
vertices (i.e., vertices connected by an edge of the graph) have different colors, with a vertex 
weighting w that either disfavors (for < w < 1) or favors (for w > 1) a given color [l|, 3]. 
Since all of the colors are, a priori, equivalent, it does not matter which color one takes to 
be given the weighting. An assignment of q colors to the vertices of a graph G, such that 
adjacent vertices have different colors, is called a "proper g-coloring" of the vertices of G. 
In the present paper we shall study a generalization of this problem in which one performs 
a proper g-coloring of the vertices of a graph G such that s colors are favored or disfavored 
relative to the remaining q — s colors. We denote these coloring problems as the DFSCP and 
FSCP for disfavored or favored weighted-set graph vertex coloring problems. We analyze 
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the properties of an associated weighted-set chromatic polynomial, denoted Ph(G,q,s,w), 
which generalizes the chromatic polynomial P(G, q) and the single-color weighted chromatic 
polynomial Ph(G,q,w) = Ph(G,q,l,w) analyzed in Ret. |. We shall denote the set of 
integers {1, s}, representing colors, as I s and the orthogonal complement {s + 1, q} as 
Ig. To each proper g-coloring of the vertices of a graph G there corresponds a term w ns , where 
n s denotes the number of vertices assigned a color in I s . The sum of such terms resulting 
from all of these proper g-colorings of the vertices of G is the function Ph(G, q, s, w). As we 
shall show below, this is a polynomial not only in w, but also in q and s. This polynomial 
constitutes a w-dependent measure, extended from the integers to the real numbers, of the 
number of proper g-colorings of the vertices of G. In the weighted-set graph coloring problem 
for a given graph G, with q G N + being the number of colors, Ph(G, q, s, w) is a map from 
(q, s, w) G N+ x I s x [0, oo) to R. One can formally extend the domain of each of the variables 
q, s, and w to R or, indeed, C, and the latter extension is necessary when one analyzes the 
zeros of Ph(G, q, s, to). The polynomial Ph(G, q, s, w) is equivalent to the partition function 
of the g-state Potts antiferromagnet on the graph G in a set of external magnetic fields, in 
the limit where the effective spin-spin exchange coupling becomes infinitely strong, so that 
the only spin configurations contributing to this partition function are those for which spins 
on adjacent vertices are different 

Qi a. 

There has been continuing interest in the Potts 
model and chromatic and Tutte polynomials for many years; reviews of the Potts model 
include and reviews of chroma* and Tntte polynomials include |-Q 

There are several motivations for this study, arising from the areas of mathematics, 
physics, and engineering. One motivation is the intrinsic mathematical interest in graph 
coloring problems and the fact that there seems to have been very little previous study of 
weighted-set graph coloring. A second one stems from the equivalence to the statistical 
mechanics of the Potts antiferromagnet in a set of magnetic fields that disfavor or favor a 
corresponding set of spin values. A third reason for interest in this subject is the fact that 
these weighted-set graph coloring problems have practical applications. For example, the 
weighted graph coloring problem with < w < 1 (i.e., the DFSCP) describes, among other 
things, the assignment of frequencies to commercial radio broadcasting stations in an area 
such that (i) adjacent stations must use different frequencies to avoid interference and (ii) 
stations prefer to avoid transmitting on a set of s specific frequencies, e.g., because these 
are used for data-taking by a nearby radio astronomy antenna. The weighted graph coloring 
problem with w > 1 (i.e., the FSCP) describes this frequency assignment process with a 
preference for a set of s frequencies, e.g., because these are most free of interference. We 
shall especially emphasize the connections with the first two of these areas in this paper. 

We note some special cases of the weighted-set chromatic polynomial. Let us consider a 
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graph G = (V,E), defined by its set of vertices V and edges (= bonds) E. We denote the 
numbers of vertices and edges of G as n(G) = n and e(G). The values s = and w = 1 
correspond to the usual unweighted proper g-coloring of the vertices of G, so Ph{G, q, s, w) 
reduces to the usual chromatic polynomial counting the number of proper g-colorings of the 
vertices of G: 

Ph(G, q, s, 1) = Ph(G, q, 0, w) = P(G, q) . (1.1) 

Since the right-hand side of Eq. (II. ip is independent of s and w, this relation also implies 
the differential equations 

dPh(G,q,s,l) 
ds 

and 

dPh(G,q,0,w) 
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aw • (L3) 

For w — 0, one is prevented from assigning any of the s disfavored colors to any of the 
vertices, so that the problem reduces to that of a proper coloring of the vertices of G with 
q — s colors, without any weighting among them. This is described by the usual (unweighted) 
chromatic polynomial P(G, q — s), so 

Ph(G,q,s,0) = P(G,q-s) . (1.4) 

Thus, the DFSCP, described by Ph{G, q, s, w) may be regarded as interpolating between 
P{G, q) and P{G, q — s) as w decreases through real values from w = 1 to w = 0. (The case 
of no weighting, w = 1, may be considered to be the border between the DFSCP and FSCP 
regimes.) If s = q, so that all of the colors receive the same weighting, then, as is clear from 
its definition, the weighted-set chromatic polynomial reduces to w n times the unweighted 
chromatic polynomial: 

Ph(G,q,q,w) = w n P(G,q) . (1.5) 

Thus, for s = and s = q, Ph(G,q, s,w) reduces to 1 and w n times P(G,q), respectively, 
while for other values of s, in particular, for integer s in the interval 1 < s < q — 1, 
Ph(G, q, s, w) is a new polynomial which is not, in general, reducible to P(G, q). Hence, while 
retaining the term DFSCP, we shall often focus on the new cases where w lies strictly between 
1 and 0. As we shall show below, the weighted-set chromatic polynomial Ph(G,q, s,w) 
satisfies a basic symmetry relation involving the interchange of s with q — s, so that a 
knowledge of the weighted-set proper g-coloring of the vertices of a graph G with a set of s 
colors is equivalent to a knowledge of the proper coloring of the vertices of G with a set of 
q — s colors. 

There are important differences between the case s — 1 studied previously in Ref. 
and the cases 2 < s < q. For s = 1, as w increases above 1 to large positive values, the 
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favored weighting of one color is increasingly in conflict with the strict constraint that no 
two adjacent vertices have the same color. Hence, this involves competing interactions and 
frustration. In contrast, in the FSCP regime with s > 2, depending on the graph G, one 
may avoid this conflict and the resultant frustration. Specific differences will be apparent 
in our explicit results. For example, in our general result for Z(G, q, s, v, w) for the circuit 
graph C n in Eq. (15. 5p below, one term vanishes identically in the case s = 1 but is present 
for other values of s in the interval I s . (By the s <-> q — s symmetry in Eqs. f!2.4[) and fl 2 . 1 3 j) . 
this also means that another term vanishes identically for s = q — 1 but is present for other 
values of s G I s .) 

II. SOME BASIC PROPERTIES 

A. Connection of Ph(G,q, s,w) with Statistical Mechanics 

It is useful to see how the function Ph(G, q, s, w) arises in a more general statistical 
mechanical context. As before, we have G = (V,E). A spanning subgraph G' C G is 
defined as the subgraph containing the same set of vertices V and a subset of the edges of 
G; G' = (V, E') with E' C E. We denote the number of connected components of G as 
k(G) and the connected subgraphs of a spanning subgraph G' as G' { , i = l,..,k(G'). To 
obtain an expression for Ph(G,q, s,w), we make use of the fact that it is a special case of 
the partition function for the g-state Potts model in the presence of external magnetic fields 
in the limit of infinitely strong antiferromagnetic spin-spin coupling. In thermal equilibrium 
at temperature T, the general Potts model partition function is given by 

{cr„} 

with the Hamiltonian 

q 

n = -jJ2 s*w - E H ? E ' ( 2 - 2 ) 

(ij) p=l i 

where i, j, i label vertices of G, Oi are classical spin variables on these vertices, taking values 
in the set I q — {1, q}, /3 — (A; B T) _1 , (ij) denote pairs of adjacent vertices, p is an integer, 
p 6 I q , and H p is an external magnetic field given by 

f H ^ for 1 < p < s 
H p = \ . (2.3) 

I for s + 1 < p < q 

The zero-field Potts model Hamiltonian 7i and partition function Z are invariant under the 

global transformation in which Oi — > gai V i G V, with g G S q , where S q is the symmetric 
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(= permutation) group on q objects. Because of this invariance, we can, without loss of 
generality, take the external magnetic fields H p to single out a set of s contiguous spin values 
(equivalently, colors) oi G I s as disfavored or favored, relative to the orthogonal complement 
of values ffj 6 if. In the presence of the magnetic fields H p given in Eq. (I2.3p . the symmetry 
group of TC and Z is reduced to the tensor product 

S q -> S s ® V* • (2-4) 

That is, if g\ G S g and g 2 G S q _ s , then the global transformation Oi — > (<?i <S> <?2)o"i V z leaves 
H and Z invariant. Here (pi ® ^2)^ means gicTj if crj G J s and if G Ijf. 
Let us introduce the notation 

K = /3J , h = (3H , y = e^, u = y - 1 . w = e' 1 . (2.5) 

Thus, the physical ranges of v are v > for the Potts ferromagnet, and — 1 < t> < for the 
Potts antiferromagnet. The weighted-set chromatic polynomial is then obtained by choosing 
the antiferro magnetic sign of the spin-spin coupling, J < and taking K — > — 00 while 
keeping h = j3H fixed. Since K = j3J, the limit K — > —00 results if (i) one takes J — > —00 
while holding T and if fixed and finite, or (ii) one takes T — > 0, i.e., /3 — > 00, with J fixed 
and finite and H —>■ so as to keep ft = f3H fixed and finite. The limit K — > — 00 guarantees 
that no two adjacent spins have the same value, or, in the coloring context, that no two 
vertices have the same color. One sees that in this statistical mechanics context, it is the 
external magnetic fields that produce the weighting that favors or disfavors a given set I s of 
spin values. Positive H gives a weighting that favors spin configurations in which spins have 
values in the set I s , or equivalently, vertex colorings with colors in this set, while negative 
H disfavors such configurations. For positive and negative H, the corresponding ranges of 
w are w > l^and < w < 1, respectively. 

In Ref. Q a formula was derived for the partition function Z which does not make any 
explicit reference to the spins Oi or the summation over spin configurations, but instead 
expresses this function as a sum of terms arising from the 2 e ^ spanning subgraphs G' C G, 
namely 

fc(G') 

Z(G, q, s,v,w) = ^2 v 6(G,) II (q~s + sw n . (2.6) 

G'CG i=l 

This generalizes a spanning subgraph formula for Z in the case s = 1 due to F. Y. Wu 



17| . which, itself, generalized the Fortuin-Kasteleyn formula for the zero- field Potts model 
partition function, 

Z(G,q,v)= ;>>^V (G,) . (2.7) 



G'CG 
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The original definition of the Potts model, ( 12. X[) and (12.21) . requires q to be in the set of 
positive integers N + and s to be a non-negative integer. These restrictions are removed by 
Eq. ( 12.61) . Furthermore, Eq. ( 12. 61) shows that Z is a polynomial in the variables g, s, v, and 
u>, hence our notation Z(G, q, s, v, w). 

The K — > — oo limit that yields the weighted-set chromatic polynomial is equivalent to 
u = —1, so 

Ph(G, q, s, w) — Z(G, q, s, — 1, w) . (2.8) 

Hence, 

fc(G') 

Ph(G,q,s,w) = ^(-l) e(G,) J] (g - s + sw n(G > } ) . (2.9) 

G'CG i=l 

We recall the factorization 

m— 1 

w m - 1 = {w- l)J2 wj ( 2 - 10 ) 

i=o 

and apply it to Eq. (12. 6p with m = n(G-). Since the variable s only appears in Eq. (12.61) in 
the form 

fc(G') fc(G') "(G-)-l 

n (g- s + sw n( ^) = n ^+ s ( w -i) j2 wr ) > ( 2 - n ) 

j=l j=l r=0 

it follows that Z(G, q, s, v, w) and Ph(G, q, s, w) can equivalently be written as polynomials 
in the variables q, v, w, and 

t = s(w-l). (2.12) 

The advantage of doing this is that it shortens expressions for these polynomials; however, 
it renders the symmetries (12.131) and (12.141) below not manifest in the resultant expressions. 

Having shown the connection of Ph(G,q, s,w) to Z(G,q, s,v,w), we observe that vari- 
ous properties of Ph(G, q, s, w) can be expressed more generally as corresponding proper- 
ties of Z(G,q, s,v,w). From Eq. (12.61) it follows that the Potts model partition function 
Z(G, q, s, v, w) satisfies a basic symmetry relating the values s and q — s: 

Z(G, q, s, v, w) = w n Z(G, q,q — s, v, w~ x ) (2-13) 

so that, in particular, setting v = —1, 

Ph(G, q, s, w) = w n Ph(G, q,q-s, w' 1 ) . (2.14) 

The symmetry relation (I2.13P is obvious from a statistical mechanics context as well as from 
the formula Eq. (12. 6p ; it is a statement of the fact that the presence of the magnetic field 
disfavors or favors the set of spin values ai G I s relative to the orthogonal complement of 



6 



spin values cr, G Ij~ , but, up to the prefactor, this is equivalent to replacing s by g — s and 
reversing the sign of H, i.e., replacing w by 1/w. 

If the magnetic field is zero, i.e., w = 1, or s = 0, so that no spin values are weighted 
differently by this field, we have 

Z(G, q, s, v, 1) = Z(G, q, 0, v, w) = Z(G, q, v) , (2.15) 

so that 

0Z(G,g,s,M) _ Q 

and 

dZ{G,q,Q,v,w) _ Q 

If s = g, so that all spin values receive the same weighting, then 

Z(G,q,q,v,w) = w n Z(G,q,v) . (2.18) 

Note that this result also follows by applying the symmetry relation (12.131) . so that 
Z(G,q,q,v,w) = w n Z(G,q,0,v,w~ 1 ) = w n Z(G,q,v). Moreover, if the disfavoring is to- 
tal, i.e., w — 0, then 

Z(G, q, s, v, 0) = Z(G, q - s, v) . (2.19) 

From the definition of Ph(G, q, s, w) as a sum of terms w na corresponding to proper 
g-colorings of the vertices of the graph G such that n s vertices are assigned colors in 
the weighted set I s , we can infer a general inequality. Let us denote the total set of 
proper g-colorings of the vertices of G as {a} and a subset as {a}subset, and let us define 
Ph(G, q, s, w) su bset as the sum of terms w ns resulting from the contributions of the subset 
{v}subset of proper g-colorings of the vertices of G. Then, since for the coloring problem 
at hand, where w > 0, each such proper g-coloring contributes a non- negative term to 
Ph(G, q, s, w), we have the general inequality 

Ph{G, q, s, w) > Ph{G, q, s, w) subset . (2.20) 

B. Some Properties Connected with Characteristics of Graphs 

If G has any loop, defined as an edge that connects a vertex to itself, then a proper 
g-coloring is impossible. This is because such a g-coloring requires that any two adjacent 
vertices have different colors, but since the vertices connected by an edge are adjacent, the 
presence of a loop in G means that a vertex is adjacent to itself. Thus, Ph(G, q, s, w) — if 
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(2.16) 
(2.17) 



G contains a loop. Hence, with no loss of generality, in our discussions of Ph(G, q, s, w) we 
shall restrict our analysis in this paper to loopless graphs G. Thus, in the text below, where 
G = (V, E) is characterized as having a non-empty edge set E ^ 0, it is understood that E 
does not contain any loops. 

Another basic property of a chromatic polynomial is that as long as two vertices are joined 
by an edge, adding more edges connecting these vertices does not change the chromatic 
polynomial. This is clear from the fact that the chromatic polynomial counts the number 
of proper g-colorings of the vertices of G, and the relevant condition - that two adjacent 
vertices must have different colors - is the same regardless of whether one or more than one 
edges join these vertices. Let us define an operation of "reduction of multiple edge(s)" in 
G, denoted Re(G), as follows: if two vertices are joined by a multiple edge, then delete all 
but one of these edges, and carry out this reduction on all edges, so that the resultant graph 
Re{G) has only single edges. Then if G is a graph that contains one or more multiple edges 
joining some set(s) of vertices, P(G,q) = P(i2^(G), 5). Since the same proper g-coloring 
condition holds for the weighted-set chromatic polynomial, we have 



Moreover, if G consists of two disjoint parts, G\ and G2, then Z(G,q, s,v,w) is simply 
the product Z(G,q, s,v,w) = Z(Gi,q, s,v,w)Z(G2,q, s,v,w), and the same factorization 
property holds for the special case v — —I that yields Ph(G,q, s,w). Hence, without loss 
of generality, unless otherwise indicated, we shall restrict our discussion here to connected 
graphs G. 



Next, we prove some general structural properties of Z(G, q, s, v , w) and Ph(G,q, s,w) 
that hold for an abitrary graph G. From Eqs. (12.61) and (12. 9p one can derive certain factor- 
ization properties of these polynomials. It is convenient to define the notation 



From Eq. (12. 6j) . it follows that we can write Z(G, q, s, v, w) in several equivalent ways: 



Ph(G, q, s, w) = PIi(Re(G), q, s, w) . 



(2.21) 



C. Properties of Coefficients in Polynomial Expansions 



q = q — s . 



(2.22) 



n 



e(G) 



n 



e(G) 




i,j,£=0 k=0 



i,j,£=0 k=0 
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n e(G) n e(G) 

= Yl c i , jtk/ q i s j v k w i = ^2 dij^qWw 1 , (2.23) 

i,j,£=0 k=0 i,j,£=0 k=0 

where ai,j,k,z, h,j,k,e, c i,j,k,i, and d^j^t are integers (and i, j, k, £ are dummy summation 
variables here), Some cn t j,k,e and can be negative, but the nonzero Ci,j,k/ and dij^j 

are positive, as follows from Eq. (12. 6p and Eq. (12.1 ip . From these equations, one infers 
corresponding ones for Ph(G, q, s, w) by setting v = —1, i.e., y — 0. 

For our analysis below and for comparisons with chromatic polynomials, three types 
of polynomial expansions will be useful. Because of the basic symmetry (12.131) . the most 
useful expansion of Z(G,q, s,v,w) is as a sum of powers of w with coefficients, denoted as 
/3a,j(q, s, v), which are polynomials in q, s, and v: 

n 

Z(G, q, s, v,w) = J2 Pz,G,j (q, s, v) w 3 . (2.24) 

3=0 

The symmetry (12.131) implies the following relation among the coefficients: 

Pz,G,j(q, s, v) = p z ,G,n-j(q, q~ s, v) for < j < n . (2.25) 
In particular, for the special case v — — 1 of primary interest here, we write 

n 

Ph(G, q, s,w) = J2 PgAQ, s ) wj > ( 2 - 26 ) 

3=0 

where 

PGAQ,s)=p ZtGjj (q,s,-l) (2.27) 

From (I2.25p . we have 

PGA9> s ) = PG*- j (q,q-s) for0<j<n. (2.28) 

If n is even, say n = 2m, then the middle coefficient is transformed into itself, giving rise to 
the results that 

If n = 2m is even, then Pz,G,m(q, s, v) = f3 z , G , m (q, q- s,v) 

P G ,m{<l, s )= fornix Q~ s ) ( 2 -29) 

From Eq. (12.61) . it is clear that the term of highest degree in w arises from products of 
the sw n ( G '^ factors over the various connected components G\ for each spanning subgraph 
G' C G, and then over the spanning subgraphs G'. This product does not involve q, so that 

Pz,G,n{ ( li s i v ) anc ^ ftG,n(q, s ) are independent of q . (2.30) 
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The /3z,G,j(q, s,v) coefficients have especially simple factorization properties, which we ana- 
lyze next. Evaluating Eq. (I2.24p at w = 0, where only the w° term remains, and combining 
this evaluation with the relation ( 12.191) . we derive the result 

(3 z , G Aq,s,v) = Z(G,q-s,v) . (2.31) 

Combining the relation for w = 1 in Eq. (12.151) with Eq. (I2.24p . we derive a formula for 
the sum of the coefficients (3z,G,j{o.i s -> v )'- 

n 

^2l3 Z)Gd (q,s,v) = Z(G,q,v) . (2.32) 

3=0 

Because this sum is independent of s, Eq. (I2.32p also yields the differential equation 

n 

^£/W3. «,«) = «• (2-33) 

3=0 

Next, we set s = q in Eq. (12.241) and use Eq. (I2.18p . Since the resulting expression must 
be proportional to w n , all of the coefficients of the terms in Z(G, q, s, v, w) of lower degree 
in w than n must vanish. Because these coefficients f3z,G,j{q, s , v ) are polynomials in q and 
s (as well as v), this means that they must contain the factor (q — s): 

Pz,G,j(q, s ) v ) an d @G,j{q, s ) contain the factor (q — s) for < j < n — 1 . (2.34) 

Furthermore, using Eq. (I2.18P for this s = q evaluation, we infer that Pz,G,n(q,q,v) = 
Z(G,q,v). However, since, by Eq. (12.301) . (3z,G,n(q, s,v) is independent of q and is only a 
function of s and v, this implies that 

(3z,G,n{q, s, v) = Z(G, s, v) (2.35) 

(so we could drop the argument q, but for uniformity with other coefficients 0z,G,j(q, s j v ), 
we shall retain it). 

Now setting s = reduces Z(G,q, s,v,w) to Z(G,q,v) (cf. Eq. (12.151) ). Since the w° 
term (given in Eq. (I2.3ip ) is, by itself, equal to Z(G,q,v) for s = 0, this means that all 
of the other terms proportional to nonzero powers , j = l,...,n in Z(G,q, s,v,w) must 
vanish when s = 0. This proves that 

For 1 < j < n, /3z,Gj{q, s, v) and Pcj(q, s ) contain a factor of s . (2.36) 

Various special cases of these results for the weighted-set chromatic polynomial are ob- 
tained by setting v — — 1 in the requisite equations. Thus, Eq. (12.311) implies 

/3 Gt0 (q,s) = P(G,q-s) , (2.37) 
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and Eq. ( 12. 35ft implies 

/3 G , n (q,s) = P(G,s) . (2.38) 

We now focus on Ph(G,q,s,w). The chromatic number of G, denoted x(G), is the 
minimal number of colors for which one can carry out a proper g-coloring of the vertices of 
G. Since the proper g-coloring constraint cannot be satisfied for integer q in the interval 
< q < x(G) — 1, the chromatic polynomial P(G,q) vanishes for these values and hence 
contains n^i?" 1 ^ ~~ 3) as a f ac t° r - Applying this to Eq. (12.381) shows that 

X(G)-1 

0G,n(q, s ) contains the factor (s — j) , (2.39) 

j=o 

and applying it to Eq. (12.371) . taking into account the shift s — > q — s, shows that 

X(G)-1 

Pco(q, s ) contains the factor j J (q — s — j) . (2.40) 

3=0 

In particular, provided that G = (V, E) contains at least one edge, so that x{G) > 2, we 
have the results 

If E ^ 0, then f3c, n (q, s) contains a factor s(s — 1) (2-41) 

and 

If E 0, then Pco(q, s) contains a factor (q — s)(q — s — 1) . (2.42) 

Thus, although the maximal degree of Ph(G, q, s, w) in w is, in general, n, it is less than n if 
s = or s = 1. In the s = case, all dependence on w disappears (cf. Eq. (11.11) ). while for 
s = 1 we have previously analyzed the maximal degree of Ph(G, q,l,w) for various families 
of graphs in Ref. For s = 1, provided that G contains at least one edge, Ph(G, q,w) 
has a factor of (q — 1) [2]. This is not, in general, true for s in the interval 2 < s < q — 1, 
and this is one of the ways that the properties of Ph(G, q, s, w) for 2 < s < q — 1 differ from 
those for s = 1. 

One can also express Z(G,q,s,v,w) as a polynomial in g with coefficients, denoted as 
az,G,i( s , v , w ), which are polynomials in s, v, and w: 

n 

Z(G, q, s, v,w) = y^ y a z ,G,n-j(s, v, w) q n ~ J . (2.43) 

3=0 

Accordingly, with the notation 

a G , n -j(s,w) = a z ,G,n-j(s, -l,w) , (2.44) 
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we write 

n 

Ph(G, q, s,w) = J2 <XG,n-j(s, w) q n ~ j . (2.45) 

j=0 

From our discussion above, we have the following results for these latter coefficients for 
Ph(G, q, s, w): 

a G , n (s,w) = l (2.46) 

and, using also Eq. f)2.2ip . 

a G , n - 1 (s,w)=ns(w-l)-e(R E (G)) . (2.47) 

Finally, since the variable s only enters Z via the combination t = s(w — 1), it is also useful 
to express the coefficients in Eq. (12.431) as polynomials in t, v, and w, and the coefficients 
in Eq. (I2.45P as polynomials in t and w. For a given graph G, we find that this usually 
simplifies the expressions. 

A chromatic polynomial P(G, q), written in the form 

n-k{G) 

P{G,q)= Yl a G,n-jq n - J , (2.48) 

3=0 

has the property that the signs of the coefficients a^n-j alternate: 

sgn(a Gin _ J ) = (-iy , < j < n - k(G) . (2.49) 

(where, as before, k(G) denotes the number of components of G, and we shall continue, with- 
out loss of generality, to focus on connected graphs, so that k(G) = 1.) This sign alternation 
property can be proved by iterated application of the deletion-contraction relation. Since 
the weighted-set chromatic polynomial Ph(G,q,s,w) does not, in general, obey a deletion- 
contraction relation, except for the values w — 1, w — 0, and s = for which it reduces to a 
chromatic polynomial, one does not expect the corresponding coefficients ac, n -j(s, w) in Eq. 
(I2.45P to have this sign- alternation property in general, and they do not. However, we have 
proved that if w is in the DFSCP interval < w < 1, then the sign alternation property 
again holds, i.e., 

sgn(a Gi „_ i (s, w)) = for < w < 1 and < j < n - 1 (2.50) 

The technical mathematical details of our proof will be given elsewhere. For the borderline 
cases w = 1 and w = 0, as well as for s = and s = q, Ph(G, q, s, w) reduces to a chromatic 
polynomial, so the sign-alternation property is already established. For j = n, namely for the 
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q° term in Ph(G, q, s, w), the sign alternation also holds for < w < 1; here the coefficient 
a G.o{ s 7 w ) contains the factor t = s(w — 1) and hence vanishes at w — 1 and s = 0. 

Setting q = in Eq. ( 12.61) . and recalling the factorization in Eq. ( 12. lip , we deduce that 

Z(G, 0, v , s, w) = az,G,o(s, v, w) contains a factor of t = s(w — 1) . (2-51) 

The same holds, a fortiori, for the v — — 1 special case, Ph(G, 0, s, w), i.e., ac,o{s, w) contains 
a factor of t = s(w — 1). 

From a study of chromatic polynomials, R. Read observed that the magnitudes of the 
coefficients of successive powers of q n ~\ < j < n — k(G) in a chromatic polynomial satisfy 
a unimodal property |8]. That is, the magnitudes of these coefficients get successively larger 
and larger, and then smaller and smaller, as j increases from to n — k(G). There is thus a 
unique maximal-magnitude coefficient, or two successive coefficients whose magnitudes are 
equal. From our calculations of weighted chromatic polynomials for a number of families of 
graphs, we have observed that in the interval < w < 1 this property continues to hold. 
We therefore state the following conjecture: Conject. Let Ph(G, q, s, w) be written as in Eq. 
(I2.45p . Then for real w in the interval < w < 1, the quantities (— l)%G,n-i(s, w), < j < 
n, are positive and satisfy the unimodal property, i.e., (— iyaG, n -j(s, w) get progressively 
larger and larger, and a maximal value is reached for a given j, or for two successive j values, 
and then the quantities (— l)i oiG,n-j{.s, w) get progressively smaller, as j increases from to 
n. 

D. Measure of Deviation from Deletion-Contraction Relation 

For a graph G, let us denote the graph obtained by deleting an edge e € E as G — e 
and the graph obtained by deleting this edge and identifying the two vertices that had been 
connected by it as G/e. The Potts model partition function satisfies the deletion-contraction 
relation (DCR) 

Z(G, q, v) = Z(G - e, q, v) + vZ(G/e, q, v) (2.52) 
and, setting v — — 1, the chromatic polynomial thus satisfies the DCR 

P(G, q, v) = P(G -e,q)- P(G/e, q) . (2.53) 

However, in general, neither Z(G,q, s,v,w) nor Ph(G,q, s,w) satisfies the respective 
deletion-contraction relation. For the special cases w = 1 and s = for which Z(G, q, s, v, w) 
and Ph(G, q, s, w) reduce to Z(G, q, v) and P(G, q), and for the special case w = for which 
Z(G,q, s,v,w) and Ph(G,q, s,w) reduce to Z{G,q — s,v) and Ph{G,q — s), respectively, 
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they do satisfy deletion-contraction relation. Hence, the deviations from such a relation for 
Z(G, q, s, v, w) and Ph(G, q, s, w) vanish in these three cases. For a given G, it is of interest 
to examine the quantities that measure the deviation from the DCR, namely 

[AZ(G,e,q,s,v,w)] DC R = Z(G,q,s,v,w)- Z(G-e,q,s,v,w)+vZ(G/e,q,s,v,w) (2.54) 

and 

[APh(G, e, q, s, w)} DCR = [AZ(G, e, q, s, -1, w)] DCR . (2.55) 

From our discussion above, it follows that [AZ(G, e, q, s, v, w)]dcr = 0, and hence also 
APh(G,e,q, s,w)]dcr = 0, for w = 1, w = 0, and s = 0; therefore, since these functions 
are polynomials in these variables, they contain a factor sw(w — 1). Moreover, since the 
condition v — is equivalent to the absence of any edges, whence Z = (q + t) n and the 
deletion-contraction relation is satisfied trivially, [AZ(G, e, q, s, v, w)]dcr also vanishes for 
v — 0. Thus, in general, 

[AZ(G, e, q, s, v, w)}dcr contains the factor svw{w — 1) , (2.56) 

and [APh(G, e, q, s, w)]dcr contains a factor of sw(w — 1). As an illustration, using our 
explicit calculations for n- vertex line graphs L n and circuit graphs C n , we find the following 
results. For the first two graphs, L 3 and C3, the deletion and contraction on any edge gives 
the same result, so we need not specify which edge is involved. We find, for any edge e, 

[AZ(L 2 , e, q, s, v, w)} DCR = svw(w - 1) , (2.57) 

[AZ(L 3 , e, q, s, w)] DC r = svw(w - 1) s(w - 1) + wv + gj , (2.58) 

and 

[AZ(C 3 , e, q, s, v, w)]dcr — svw(w — 1) wv 2 + 2wv + s(w — 1) + gj . (2.59) 

As before, the corresponding [APh(G, e, q, s,w)] DC r expressions are obtained by set- 
ting v = —1 in these equations. It is straightforward to calculate similar differences 
[AZ(G,e,q, s,v,w)]dcr f° r graphs with more vertices and edges, but these are sufficient 
for our illustration. 



E. Distinguishing Between Various Equivalence Classes of Graphs 

An important property of the weighted-set chromatic polynomial Ph(G, q, s, w) is the 
fact that it can distinguish between certain graphs that yield the same chromatic polyno- 
mial P(G,q). This is true for all w and s values except the special values w — 1, w — 0, 
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s = 0, and s = q, for which Ph(G, q, s, w) is reducible to a chromatic polynomial. More 
generally, an important property of the partition function of the Potts model in a set of 
nonzero external magnetic fields of the form ( 12. 31) . Z(G,q, s,v,w), is that this function can 
distinguish between graphs that yield the same zero-field Potts model partition function, 
Z(G, q, s, v, 1) = Z(G, q, v). Two graphs G and H are defined as (i) Tutte-equivalent if they 
have the same Tutte polynomial, or equivalently, the same zero-field Potts model partition 
function, Z(G,q,v), and (ii) chromatically equivalent if they have the same chromatic poly- 
nomial, P(G, q). Here we recall that the Tutte polynomial T(G, x,y) of a graph G is defined 
as 

T(G, x,y)= J2( x ~ l) k{G ' ] ' HG \y - 1) C(G ° (2.60) 

G'CG 

where G' is a spanning subgraph of G and c{G') denotes the number of (linearly independent) 
cycles in G' . This polynomial is equivalent to the zero-field Potts model partition function, 
via the relation 

Z(G, q, v) = (x- l) k ^(y - l) n T(G, x, y) . (2.61) 
where y = v + 1 as in Eq. (I2.5P and 

x = 1 + - . (2.62) 

v 

We give some examples. Recall the definition that a tree graph is a connected graph that 
contains no circuits (cycles). The set of tree graphs with n vertices, generically denoted 
{T n }, forms a Tutte equivalence class, with T(T n ,x,y) = x 71 " 1 , or equivalently, 

Z(T n ,q,v) =q(q + v) n - 1 . (2.63) 

However, Z(G,q,s,v,w) is able to distinguish between different tree graphs in a Tutte- 
equivalence class. Because Z(G, q, s, v, w) reduces to a zero-field Potts model partition func- 
tion for w — 1, w — 0, s — 0, and s — q, it follows that the difference between Z(G, q, s, v, w) 
and Z(H,q, s,v,w) for two Tutte-equivalent graph G and H must vanish if w — 1, w — 0, 
s = 0, or s = q. This difference also vanishes for v = 0, because in this case the only 
spanning subgraph that contributes to the sum in Eq. (12.61) is the one with no edges, which 
is the same for any (connected) Tutte-equivalent G and H. Since these are all polynomials, 
it thus follows that 

Z(G, q, s, v, w) — Z(H, q, s, v, w) contains the factor s(q — s)vw(w — 1) . (2.64) 
As an illustration, using our results for Z(L±, q, s, v, w) and Z{S±, q, s, v, w), we have 

Z(Si, q, s, v, w) — Z(Li, q, s, v, w) = s(q — s)v 2 w(w — l) 2 (2.65) 
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and consequently 

Ph^S^, q, s, w) — Ph(Li, q, s, w) = s(q — s)w(w — l) 2 . (2.66) 

The zero-field Potts model partition function Z(G,q,v), or equivalently, the Tutte poly- 
nomial T(G,x,y), encodes information on the number of (linearly independent) cycles con- 
tained in the graph G, as is evident from the definition (12.601) . Define two scaled variables 
as 

q'= q -, v' = ~. (2.67) 

s s 

Let us us consider a graph, denoted G nc , which contains no cycles (nc), i.e., which has 
c(G) = 0. A connected graph of this type is a tree graph, while a general graph is called a 
forest. For a graph G nc we find the following scaling relation: 

Z(G nc , q, s, v, w) = s n Z(G nc , q' , 1, v', w) . (2.68) 

This is proved as follows. We start with the cluster formula (12 .6p and rewrite this as 

KG') 

Z(G nc , q, s, v,w)= K) e(G ° s< G '^ J] (q> - 1 + w n ^) . (2.69) 

We next use the relation, which holds for any graph G', 

c(G') + ra(G') = e(G') + fc(G') (2.70) 

and the fact that n(G') = n(G) = n to rewrite the factor s e ( G ')+ k (G') as s c{G')+n _ gj nce Q nc 
has no cycles, it follows that the same is true for any subgraph of G nc , in particular, the 
spanning subgraph G', so c(G') = 0. Hence, we can move the factor of s n in front of the 
summation, and we have 

KG') 

Z(G nc ,q,s,v,w) = s n ^') e{GI) II (q'~l + w n W) 

G (ZG nc i=l 

= s n Z{G nc ,q\l,v\w) . (2.71) 

□ 

Hence, the difference 

[AZ(G, q, s, v, w)] cydes = Z(G, q, s, v, w) - s n Z(G, q', 1, v', w) (2.72) 

provides a measure of the number of cycles in G. For example, using our general result for 
Z(C n , q, s, v, w) given below in Eq. (15.21) . we calculate 

W7tn M (s - l)(q - s + sw n )v n 

[AZ(C n , q, s, v, w)\ cydes = (2.73) 
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Moreover, the proper g-coloring condition implies that if two different graphs G and H 
differ only in having different numbers of edges connecting a pair of adjacent vertices, for 
one or more such pairs, so that Re{G) = Re(H), then Ph(G,q,s,w) = Ph(H,q,s,w). A 
simple example is provided by the line and circuit graphs with n = 2 vertices, L 2 and C 2 , 
the latter of which has a double edge connecting the two vertices. Using our results in Eqs. 
( 13. 2p and (15. 7p . we calculate the difference 

Z(C 2 , q, s, v, w) - Z(L 2 , q, s, v, w) = v(v + l)[q+ s(s - l)(w + 1)] . (2.74) 

The fact that the difference in Eq. (I2.74p vanishes for v = —1, i.e., that Ph(L 2 , q, s,w) = 
Ph(C 2 , q, s, w), is a special case of the general result (I2.2ip . 

In the context of graph coloring, since s 6 I s , if one sets q to a particular value, this 
implicitly sets a corresponding upper bound on s. In particular, if q = 1, then s can take on 
only the values or 1, and these are related by the symmetry (I2.13p . For s = 0, we have, 
by Eq. f[2TT5]) . that 

Z(G, l,0,v,w) = Z(G, l,v) = y e{ - G) w n , (2.75) 
where y = v + 1. For s = 1, applying Eq. (I2.13p . we have 

Z(G, 1, 1, v, w) = w n Z{G, 1, 0, v, w' 1 ) = w n Z{G, 1, v) = y e{G) w 2n . (2.76) 

If G has at least one edge, then the right-hand sides of both Eq. (I2.75P and Eq. (12.761) 
vanish for the case y = (v = —1) that yields the weighted-set chromatic polynomal. In 
order for two graphs G and H to be chromatically equivalent, a necessary condition is that 
they must have the same number of vertices, n(G) = n(H), since the degree in q of P(G, q) 
is n(G). An elementary property of the chromatic polynomial P(G,q), proved by iterative 
application of the deletion-contraction theorem, is that the coefficient of the g™^) -1 term is 
—c(Re{G)). Therefore, another necessary condition that two graphs G and H must satisfy 
in order to be chromatically equivalent is that e(-Rs(G)) = e(RE(H)). If G contains at 
least one edge, then Ph(G,l, s,w) = 0. Note here that since s is bounded above by q, it 
follows that if q = 1, then s can only take on the values s — or s — 1. Hence, if G 
and H are chromatically equivalent, then either (i) neither contains any edges, in which 
case Ph(G,q, s,w) = Ph(H,q, s,w) = (q + t) n , where n = n{G) = n(H), or (ii) if G, and 
hence H, contains at least one edge, then Ph(G,l, s,w) = Ph(H,l, s,w) = 0. Hence, if 
G and H are chromatically equivalent and contain at least one edge, then the difference 
Ph(G,q, s,w) — Ph(H,q,s,w) contains a factor that vanishes when q = 1. Because of the 
implicit condition on s for a given q, this factor is not, in general, (q — 1). As an example, 
the difference Ph(S^ q, s, w) — P/i(L 4 , q, s, w) in Eq. (I2.66P contains the factor s(q — s). For 



17 



q = 1, the values of s are implicitly restricted to s = and s = 1. For either of these choices, 
the factor, and hence the difference, vanishes. 

A remark concerning duality for planar graphs is also in order here. Let G = (V, E) 
be a planar graph, and denote its planar dual by G*. The chromatic polynomial P(G,q) 
counts not just the proper g-colorings of the vertices of G, but also, and equivalently, the 
proper q colorings of the faces of G*. Similarly, for this planar graph G, the weighted-set 
chromatic polynomial Ph(G, q, s, w) describes not just the weighted-set proper g-colorings of 
the vertices of G but also, and equivalently, the weighted-set proper g-colorings of the faces 
of G*. 

F. Lower Bounds on Ph(G,q, s,w) 

We derive some bounds on Ph(G,q, s,w) for certain types of graphs. Our method for 
this will be to calculate the contribution to Ph(G,q, s,w) resulting from a certain proce- 
dure for performing proper g-colorings of the graph G. By the general formula expressing 
Ph(G, q, s, w) as the v — — 1 special case of the Potts model partition function Z(G, q, s, v, w) 
together with the formulation of this partition function as a sum over spin (or equivalently, 
color) configurations, Eqs. (12. lj) with (12.21) . it follows that there are other color configu- 
rations in addition to the particular one that we consider, contributing (positive terms) to 
Ph(G, q, s, w). Therefore each particular proper g-coloring procedure that we consider yields 
a lower bound on Ph(G,q, s,w). The specific proper g-coloring procedure that provides a 
good lower bound to Ph(G, q, s, w) depends on the type of graph G, the values of q, s, and 
w. 

Let us consider a bipartite graph G^p, defined as a graph whose vertex set V can be 
partitioned into subsets V\ and Vj> such that a vertex in V\ has edges that connect it only to 
a vertex or vertices in Vi and vice versa. An equivalent condition for a graph to be bipartite 
is that its chromatic number xi^Hp) — 2. As above, we denote the number of vertices in G, 
as n(G) = n and, further, the number of vertices in V\ and Vi as n\ and ri2, respectively. 
With no loss of generality, we label these subsets of vertices so that n\ < n^. These numbers 
rti and ri2 may be comparable or may be quite different. For example for a lattice graph such 
as the circuit graph, the square, honeycomb, simple cubic, or body-centered cubic lattices, 
with periodic boundary conditions that preserve the bipartite nature of the lattices, n\ = rii- 
However, for the star graph S n , V\ consists of the central vertex, so that n\ = 1, while Vi 
is comprised of all of the vertices on the ends of the edges forming the rays of the star, so 
ri2 = n — 1. For an S n graph with n » 1, it follows that ri2 » ri\. 

For the following, we assume that q > 2 so that a proper g-coloring of the bipartite 
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graph Gbip is possible. If w = 1, then a well-known elementary lower bound on P(Gbi P , q) 
is obtained by (i) assigning a single color to all of the vertices in V\ and (ii) independently 
choosing a color out of the remaining q — 1 for each of the vertices in V2. There are (7(9 — l) n2 
ways of doing this. Since, in general, there are also other color configurations contributing 
to a proper g-coloring of Gbi P , this yields the lower bound 

P(G bip ,q)>q(q-l) n2 . (2.77) 

For the weighted-set chromatic polynomial, the situation is more complicated. With no 
loss of generality, we again label the vertex subsets so that n\ < n 2 . We also take q > 2 so 
that a proper (weighted) g-coloring is possible and also assume that s > 2 (and s < q, as 
discussed above). Then for sufficiently large w > 1 in the FSCP interval, one lower bound 
on Ph(Gu P , q, s, w) is obtained by (i) assigning one color from the favorably weighted set I s 
to all of the vertices in V\, and then (ii) independently, for each vertex in V2, assigning a 
color from among the remaining s — 1 colors in I s . This combined color assignment can be 
made in s(s — l)™ 2 ways, and yields a contribution s(s — l) n2 w n to Ph(Gbi P ,q, s,w). From 
the inequality (12.201) . it then follows that 

Ph{G Up , q, a, w) > s(s - l) n2 w n . (2.78) 

(Since n\ < ri2, this is an equivalent or better bound than the one obtained by making the 
above color assignments with V\ and V2 reversed, viz., Ph(Gu p ,q, s,w) > s(s — l) ni w n .) 

On the other hand, for w in the DFSCP interval < w < 1, it can be preferable to 
minimize the number of vertices with colors in I s in order to maximize the contribution to 
Ph(Gbi P , q, s, w). For sufficiently small (positive) w, provided that q > s + 2, a lower bound 
on Ph(Gbi P , q, s, w) is then obtained by (i) assigning a single color from to all of the 
vertices of V±, and (ii) independently, for each vertex of V2, assigning a color from among 
the remaining (q — s — 1) colors in ij-. This combined color assignment can be made in 
(q — s)(q — s — l)™ 2 ways. Invoking the inequality (12.201) again, we have 

Ph(G bip , q, s,w)>(q-s)(q-s- l)™ 2 . (2.79) 

If w is only slightly less than unity, say 1 — e < w < 1 for sufficiently small positive e, 
provided also that q > s + 1, a different type of lower bound on Ph(Gbi P ,q, s,w) can be 
obtained by the following proper g-coloring procedure: (i) one assigns a single color from I s 
to all of the vertices of V\ and (ii) independently for each vertex in V2, one assigns a color 
from 1^-. There are s(q — s) n2 ways of doing this, and the resultant term in Ph(Gbi P , q, s, w) 
is sw ni (q — s) 112 . Using the inequality f)2.20p again, one thus infers the lower bound 

Ph(G Mp , q, s, w) > sw ni (q - s) n2 . (2.80) 
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Which of these lower bounds is the best depends in detail on Gu p (in particular, on n\ and 
n 2 ), q, s, and w. It is easy to generalize these lower bounds to multipartite graphs. 

III. CALCULATIONS OF Z(G,q,s,v,w) AND Ph(G,q,s,w) FOR SOME FAMILIES 

OF GRAPHS 

In this section we give some illustrative explicit calculations of Z(G,q,s,v,w) and 
Ph(G, q, s, w) for various families of graphs. Although we generally consider connected 
graphs, we note that for the graph N n consisting of n vertices with no edges, 

Z{N n , q, s, v, w) = Ph(N n , q, s, w) = (q + t) n . (3.1) 



A. Line Graphs L n 

The line graph (also called path graph) L n is the graph consisting of n vertices with each 
vertex connected to the next one by one edge. In general, az,L n ,n-i(q, s, v ) — nt + (n — l)v. 
We proceed to give some explicit results for Z(L n ,q,s,w) for various values of n. The case 
Li = Ni is already covered by Eq. (13. ip . For the first few n values, we also give the 
expansions in terms of powers of q and, for this latter expansion, we use the variables q, t, 
and w instead of q, s, and w, because this makes the expressions shorter: 

Z(L 2 , q, s, v, w) = s(s + v)w 2 + 2s(q — s)w + (q — s)(q — s + v) 



q 2 + (2t + v)q + tt + v(w + 1) 



(3.2) 



Z(L°3, q, s,v,w) = s(s + v) w + s(q — s)(3s + 2v)w 



+ s(q — s) 3(q — s) + 2v w + (q — s)(q — s + v 



q 3 + (3t + 2v)q 2 + (3t 2 + 2vtw + Avt + v 2 )q 



+ t(v 2 w 2 + 2vtw + wv 2 + t 2 + 2vt + v 2 ) 



(3.3) 



For L 4 we give only the expansion in powers of w, since the equivalent expansion in powers 
of q becomes somewhat lengthy: 



Z{L&, q, s, v, w) = s(s + v) 3 w 4 + 2s(q - s)(s + v)(2s + v)w' d 
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+ s(q - s) - 3(s 2 + (q- s) 2 ) + 3q(q + v) + 2v 2 



w 



+ 2s{q - s)(q - s + v)[2(q - s) + v]w + (q - s)(q - s + v) 3 . (3.4) 



B. Star Graphs S n 

A star graph S n consists of one central vertex with degree n — 1 connected by edges 
with n — 1 outer vertices, each of which has degree 1. (The context will always make clear 
the difference between this symbol for the n-vertex star graph and the symbol S n for the 
symmetric group on n objects.) The graph S2 is degenerate in the sense that it has no central 
vertex but instead coincides with L 2 . The graph S 3 is nondegenerate, and coincides with L 3 , 
while the S n for n > 4 are distinct graphs not coinciding with those of other families. For 
n > 2, the chromatic number is x{S n ) = 2. By the use of combinatoric coloring methods, 
we have derived the following general formula for Z(S n , q, s, v, w): 



n-l / _ 1 \ 

Z{S n , q, s, v,w) = J2 ■ ) vj (? + sw ' +1 ) (? + sw) n - l - j 



(3.5) 



where q = q—s, as given in Eq. fl2.22p . Evaluating Eq. (13.51) for v — — 1 yields Ph(S n , q, s, w). 
As an explicit example, for the graph 54, we calculate 

Z(Si, q, s,v,w) = Z(T±, q,v)w 4 + s(q — s)(4s 2 + 6sv + 3v 2 )w 3 

+ 3s(q — s)[2s(q — s) + qv]w 2 + s(q — s) 4(g — s) 2 + 6(q — s)v + 3v 2 w 

+ Z(T 4 ,q-s,v) (3.6) 
where Z(T n , q, v) was given in Eq. f)2.63p . 



C. Complete Graphs K n 

The complete graph K n is the graph with n vertices such that each vertex is connected 
to every other vertex by one edge. The chromatic number is x{K n ) — n and the number of 
edges is e(K n ) = (")• Let us introduce the compact notation xg = x9(x), where 9{x) is the 
step function from R — > {0, 1} defined as 9(x) — 1 if x > and 9{x) = if x < 0. We have 
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derived the following theorem giving a general formula for Ph(K n , q, s, w): 

n 

Ph(K n ,q,s,w) = ^ (3 Kn/ (q, s) w £ (3.7) 

where 

, \ (f-l)e (n-<-l) e 

^=(") [ II( a -^][ II (l-s-m)]. (3.8) 

^ ' j=0 m=0 

Proof: This result is proved by a combinatoric coloring argument. Accordingly, we take 
q to be a non-negative integer. The resultant Eqs. (13.71) and (13.81) allow the extension of q 
to M. (and C). First, if q < x{K n ) — n i then Ph(K n ,q, s,w) vanishes identically. Hence, we 
shall formally take q > n to begin with; once we have obtained the results (13. 71) and (13.8j) . 
it will be seen that they allow an extension of q away from this range. If s > n, then one 
can assigning n different colors to the n vertices of K n from the set I s , and this gives rise 
to a term with degree n in w. To determine the coefficient of this term, we enumerate the 
number of ways this color assignment can be made. We pick a given vertex and assign some 
color from I s to this vertex, which we can do in any of s ways. Then we go on to the next 
vertex and assign one of the remaining s — 1 colors in I s to that vertex, and so on for the 
n vertices. The number of ways of making this color assignment, i.e., the coefficient of the 
term in Ph(K n , q, s, w) of maximal degree in w, viz., to", is therefore 

n-1 
3=0 

The fact that this coefficient is P(K n , s) agrees with the v = —1 special case of the general 
result of Eq. (12.351) . Similarly, the term of order w° is obtained by assigning n different 
colors to the n vertices of K n from the orthogonal set S . By reasoning analogous to that 
given above, it follows that the number of ways of doing this is given by replacing shy q — s 
in Eq. (EES) , so 

n-1 

f3 Kn ,o(q, s) = H(q - s - j) = P(K n , q - s) . (3.10) 

3=0 

Having illustrated the logic on these two extremal terms, let us next consider the general w e 
term with < £ < n. This term arises from color assignments in which we pick I different 
colors from the set I s and assign them to I of the n vertices of K n , and then n — i different 
colors from the orthogonal complement set S 1-1 , which are assigned to the remaining n — i 
vertices. The number of ways of doing this is 

(e-i) e (n-e-i) g 
/W = [ II ( s " *)] [ II (q-s-m)]. (3.11) 

j=0 m=0 
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This proves the result in Eqs. ( 13. 7ft and 1 13. 8ft . □ 

Evidently, with the polynomial Ph(K n , q, s, w) as specified in these equations, one can 
extend q and s away from non-negative integer values. Our result in Eqs. (]3.7[) and ( 13.81) 
generalizes the result for the case s = 1 given in [2|. As is evident, for w = 1 or s = 0, 
Ph(K n , q, s, w) reduces to the (usual, unweighted) chromatic polynomial 

n-l 

P(K n ,q) = H(q-j) . (3.12) 

j=0 

A corollary of Eqs. ( 13.71) and ( 13. 8ft is that 

If s < n , then ^K n ,j(Q, s) = for s < j < n (3.13) 

and hence 

deg w (Ph(K n , q, s, w)) = min(n, s) . (3.14) 

Having calculated Ph(K n ,q,s,w), it is appropriate to discuss here another aspect in 
which the weighted-set chromatic polynomial differs from the (usual, unweighted) chromatic 
polynomial. Let us consider a graph G that has the property of being composed of the union 
of two subgraphs, G = G\ U G2, such that G\ fl G2 = K m for some m. In the rest of this 
paragraph, we assume that G has this property. Then P(G, q) satisfies the relation 

fl^ ^y . (3.15) 
P(K m , q) 

(This is sometimes called the complete-graph intersection theorem (KIT) for chro- 
matic polynomials.) In contrast, in general, Ph(G,q, s,w) is not equal to 
Ph(G\,q,s,w)Ph(G2,q,s,w)/Ph(K m ,q,s,w). This equality holds only for the four val- 
ues w — 1, w = 0, s — 0, and s = q where Ph(G, q, s, w) reduces to a chromatic polynomial. 
As a measure of the deviation from equality, we define 

IAP ,m vi - Tjurn ^ Ph(Gx,q,s,w)Ph(G 2 ,q,s,w) 

[APh(G, q, s, w)\kit = Ph(G } q, s, w) — — r . (3.16) 

1 tlyJ\ m , q, s, w J 

The vanishing of [APh(G, q, s, w )}kit for w — 1, w = 0, and s = is obvious. To show that 
this vanishes for s = q, we use the relation ( 12. 18ft and obtain 

[APh(G, q, q, w)] KIT = w n \P(G, q) - ^^/^ q) ] = . (3.17) 

Combining these results with the property that [APh(G,q, s,w)]kit is a rational function 
in its arguments, we have thus shown that 

[APh(G, q, s, w)]kit contains the factor s(q — s)w(w — 1) . (3.18) 
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We give two illustrations. The line graph L 3 has the property of being comprised of two 
L 2 graphs intersecting on L\ = K\. Using (the v = — 1 special cases of) our results in Eqs. 
(EH), (|H2}, and fl^ijl . we calculate 

[APft(L 3 , 5, »)W = " 3> r (t " T, 1)2 ■ (3-19) 

q + s(w — 1) 

Similarly, the graph L 4 can be decomposed into L 3 and L 2 subgraphs that intersect on an 
L\ = K\ graph. Using our results in Eq. (13. ip and (I3.2p - (I3.4I) . we calculate 

s(q — s)w(w — l) 2 q + s(w — 1) — (w + 1) 
[APh(L 4} q, s, w)\ KIT = -i- i . (3.20) 



g + s(iy — 1) 

A slightly more complicated case is the 4-vertex graph consisting of a box with one 
diagonal edge added. This graph has the structure of two C 3 = K 3 subgraphs intersecting 
on the diagonal edge graph, L 2 = K^- For this graph we calculate 

[APh(C 4d , q, s, w)) KIT = 2s(q - s)w(w - l) 2 [l - 1)W 1 . (3.21) 

L Ph{K 2 ,q,s,w)) 

We see that each of these differences [APh(G, q, s, w)]kit satisfies the general factorization 
property of Eq. (13. 18j) . 



IV. Z(G,q,s,v) AND Ph(G,q,s,w) FOR CYCLIC STRIP GRAPHS 

A. General Structure 

Refs. [l|, have given a general structural formula for Z(G S , L y x m, BC, q, s, v, w) on strip 
graphs G s of width L y vertices and length L x , with cyclic or Mobius boundary conditions 
(BC's) for the case s — 1. Here we discuss the generalizations to arbitrary (integer) s in the 
interval < s < q. For cyclic strip graphs G s of this type we have 

L y nzh(L y ,d,s) 

Z(G s ,L y x m,cyc.,q,s,v,w) = c (d) ^ [*z,G.,Ly,dj(.Q, s, v, w)] m , (4.1) 

d=0 j=l 

where m = L x for strips of the square and triangular lattices and m = L x /2 for strips of the 
honeycomb lattice. The coefficients are given by 

c^=c^(q)^J2(-iy( 2d T j )q d - j (4.2) 
3=0 \ 3 / 
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where q = q — s, as specified in Eq. (12.221) . The first few of these coefficients are = 1, 
gW — q — 1 = q — s — 1, = g 2 — 3q + l,etc. For Mobius strips, the switching of certain 
c (d), s, as specified for s — 1 in general in |3j, [l9j, generalizes to arbitrary s in the interval 
< s < q. Further, from Eq. (14.81) . it follows that there is only one term with d = L y , and 
we find (dropping the j subscript) 

^z,L v ,L y (q,s,v } w) = v . (4.3) 
For Ph(G s , L y x m, eye, q, s, w) = Z(G, L y x m, eye, q, s, —1, w), we have 

Ly nph(L v ,d,s) 

Ph(G, L y x m, eye, q, s, w) = c {d) ^ [^Ph,G s ,L y ,d,A^ s > w )] m ( 4 - 4 ) 

d=0 j=l 

As with Z, there is only one term with d = L y , and Eq. H4.3|) shows that this is 

^Ph,G s ,L y ,L y {q,s,w) = -1 . (4.5) 
The nzh(Ly, d, s) satisfy the identity 



J2c {d) n Zh (Ly,d,s) = q L y (4.6) 

d=0 

while the nph(L y , d, s) satisfy the identity 

Ly 

J2 ~c^n Ph {L y , d) = P[L y , q) = q{q - l)**" 1 . (4.7) 

d=0 

The reason why these identities hold for general s with the same right-hand side as for 
s = and s — 1 is that the basic coloring constraints remain the same; the only thing 
that is different for nonzero s is the weighting factors. One method of calculating the 
nzh{L y ,d, s) and nph(L y ,d, s) is to differentiate these respective equations L y times. One 
thereby obtains two respective sets of L y + 1 linear equations in the L y + 1 unknowns 
nzh{L y , d, s) and nph(L y , d, s) for d = 0,1, L y . Solving these equations determines these 
numbers n Z h(L y , d, s) and nph(L y , d, s). 

We have used the method above to calculate the nzh(L y , d, s) and riph(L y , d, s). For the 
nzh{L y , d, s) we find 

n Zh {L y , L y , s) = 1 (4.8) 
n Zh (L y , L y - 1, s) = (s + l)L y + {L y - 1) (4.9) 
n Z h{L y + 1, 0, s) = (s + l)n Zh (L y , 0, s) + n Z h{L y , 1, s) (4.10) 
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and, for 1 < d < L y + 1, 

n Zh (L y + 1, d, s) = n Zh (L y + 1, d - 1, s) + (s + 2)n Zh (L y , d, s) + n Zh (L y , d+ 1, s) . (4.11) 

Some additional results, besides the general formulas for n Z h(L y , d, s) with d = L y and 
d — L y — 1 given in Eqs. (14. 8 p and (14. 9p . are 

L y = 2: n Zh (2,0,s) = s 2 + 2s + 2 (4.12) 

L y = 3 : n Zh (3, 0, s) = s 3 + 3s 2 + 6s + 5, 

n Z h(3,l,s) =3(s 2 + 3s + 3) (4.13) 

L y = 4 : n Zh (4, 0, s) = s 4 + 4s 3 + 12s 2 + 20s + 14 , 
n Zh (4, 1, s) = 4s 3 + 18s 2 + 36s + 28 , 

n Zh (4,2,s) = 6s 2 + 20s + 20 . (4.14) 

The numbers n Z h(L y , d, s) of \ z ,L y ,d,j(q, s, v, w)'s corresponding to each cfw in the general 
Potts model partition function are reduced for the special case v — — 1 that yields the 
weighted-set chromatic polynomial. By coloring combinatoric arguments similar to those 
used in 19] and 0], we determine the n Ph (L y , d, s) as follows. The numbers n Ph (L y , d, s) are 
identically zero for d > L y , and 

n Ph (L y ,L y ,s) = 1 (4.15) 

Uph(Ly, L y - 1, s) = (s + l)Ly (4.16) 

n Ph (L y + 1, 0, s) = s n Ph (L y , 0, s) + n Ph (L y , 1, s) (4.17) 
and, for 1 < d < L y + 1, 

n Ph (L y + 1, d, s) = n Ph (L y + 1, d - 1, s) + (s + l)n Ph (L y , d, s) + n Ph (L y , d+l,s) . (4.18) 

Some additional results, besides the general formulas for n P h(L y ,d, s) for d = L y and <i = 
— 1 given in Eqs. (14.151) and (I4.16p . are 

^ = 2: 7i Pfc (2,0,s) = s 2 + s + l (4.19) 



L y = 3: n Ph (3,0,s) = s 3 + s 2 + 3s + 2 
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n Ph (3, l,s) = 3s 2 + 5s + 4 (4.20) 

L y = 4: n Ph (A,0,s) = s 4 + s 3 + 6s 2 + 7s + 4 , 
n Ph (4, 1, s) = 4s 3 + 9s 2 + 15s + 9 , 

n Ph (A,2,s) = 6s 2 + lis + 8 . (4.21) 

As we have noted above, for s = 0, Z(G, q, s, v, w) reduces to Z(G, q, v), so we focus on the 
nonzero (integer) s values in the set I s . For these values we find 

n Ph (L y , d, s) = n Zh (L y , d,s-l) + n Zh (L y - 1, d, s - 1) . (4.22) 

From our determination of the numbers nzh{L y , d, s) and n P h(L y , d, s), we next calculate 
the total numbers 

Ly 

N Zh ,L y ,s = n zh(L y , d, s) (4.23) 

d=0 

and 



N Ph>Ly , s = n Ph{L y , d, s) . (4.24) 

d=0 

From Eq. (14.221) it follows that for nonzero (integer) s values in the set I s , 

N Ph>Ly , s = N Zh ,L y ,3-i + N Zh ,L y -i,s-i • (4.25) 
Using our results in Eqs. (I4.8p - fl4.lip . we find that, for nonzero s£/ s , 



3=0 

and hence, using Eq. f)4.25p . we have 

L 



— - E ft) (?) (. - V- + E ( £ »: l ) (?) (. - u-" ■ (-) 

A few explicit results for low values of s are 

W*M,. = s + 2 ( 4 - 28 ) 
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and 



and 



N Zh ,2,s = s 2 + 4s + 6 (4.29) 

N Zh ,3,s = s 3 + 6s 2 + 18s + 20 = (s + 2)(s 2 + 4s + 10) (4.30) 

Nzha* = si + 8s3 + 36s2 + 80s + 70 ( 4 - 31 ) 

N Z h,5,s = s 5 + 10s 4 + 60s 3 + 200s 2 + 350s + 252 

= (s + 2)(s 4 + 8s 3 + 44s 2 + 112s + 126) (4.32) 

N Zh ,6,s = s 6 + 12s 5 + 90s 4 + 400s 3 + 1050s 2 + 1512s + 924 (4.33) 

N Ph ,x, s = s + 2 (4.34) 

N Ph ,2,s = s 2 + 3s + 4 (4.35) 

jV PM)S = s 3 + 4s 2 + lls + 10 (4.36) 

Np hAs = s 4 + 5s 3 + 21s 2 + 37s + 26 (4.37) 

N PhAs = s 5 + 6s 4 + 34s 3 + 88s 2 + 123s + 70 (4.38) 

Np hfitS = s 6 + 7s 5 + 50s 4 + 170s 3 + 366s 2 + 401s + 192 . (4.39) 



For L y » 1, these total numbers have the following dominant asymptotic exponential 
growth rates (suppressing power-law prefactors): 

N Zh ,L y ,s = {s + A) 1 * forL, ^ oo (4.40) 

and 

N Ph ,L y ,s = (s + 3)^ forL.^oo. (4.41) 

We note that the s — > q — s symmetry (I2.13P is not manifestly evident in the various results 
that we have given for n Zh (L y ,d, s), n Ph (L y , d, s), N Zh ,L y ,s, and N Ph)Ly>s . This symmetry 
arises via identities involving the cS d > and the A's. The symmetry (12.131) implies that (for 
s G I s and denoting L x = m) 

L y n zh (L y ,d,s) 

d=0 j=l 

Ly n zh (L y ,d,q-s) 

w n J^c^ d \s) J2 [^AifaQ-MW- 1 )]™ (4.42) 

d=0 j=l 
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and, for v = — 1, 

L y n Ph (L y ,d,s) 
L B n Ph (L y ,d,q-s) 

w n Y c(d) ( s ) E [Ap^^^g-^^- 1 )]" 1 . (4.43) 

d=0 j=l 

In particular, when s = q, Z simplifies to a multiple times the zero-field Z as speci- 
fied by the general relation (12.181) . with a consequent reduction in the number of distinct 
^z,L y ,d,j(q, s, v, w)'s. Using the fact that c^(0) = (— l) d [191], we can express this reduction 

as 

L y n zh (L y ,d,q) 

E ^Z,L yA M^^,w)} m = 
d=0 j=l 

L y n z (L y ,d) 

w n J2cW(q) Y I^zmMvT (4-44) 

d=0 j=l 

where nz(L y , d) = nzh(L y , d, 0). For s — q — 1, Eq. (I4.42[) yields the identity 

L y n zh (L y ,d,q-l) 

E c(d) d) E [Az^ |( y(?,g-l,«,«;)r = 

d=0 j=l 

L y n zh (L y ,d,l) 

w^^^-l) E l^MJ&h^w- 1 )]™ (4.45) 

d=0 j=l 



where c^(l) takes on the values 19] 

{1 if d = mod 3 
if d = 1 mod 3 . (4.46) 
-1 if d = 2 mod 3 

Analogous identities follow from Eq. (I4.42[) for s = q — 2, s = q — 3 and s = g — 4 
(assuming that s E 7°), where the values of c^(2), c^(3), and c^(4) were given in Eqs. 
(2.19)-(2.21) of Ref. [19J. Thus, in using the results above for nzh{L yi d, s), nph{L yi d, s), 
Nzh,L y ,s, an d N Ph)Ly)S1 it is understood that they apply for generic values of s G J s but 
involve simplifications for special values of s. We shall give an example of this in the next 
section. 
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V. CIRCUIT GRAPHS C n 

The circuit graph C n , or equivalently, the ID lattice with periodic boundary conditions, 
has chromatic number 

{2 if n > 2 is even 
(5.1) 
3 if n > 3 is odd 

(The case n = 1 is a single vertex with a loop, for which there is no proper g-coloring, so 
Ph(Ci, q, w) vanishes identically.) In general, Z(C n , q, s, v, w) has the structure 

s+l 

Z(C n , q, s, v, w) = J> Z)1|0J (<Z, s, v, w)] m + ~c^v m (5.2) 

where we recall that We find that (suppressing arguments) the \z,ifl,j for 

1 < i < s + 1 are given by 



\ _ 1 

Az,i,o,i - 2 



1/2 

q — s + v + w(s + v) ± {g — s + f + w(s + f )} 2 — 4w?(g + t>) 



(5.3) 



for j = 1, 2 (which is the total set for s = or s = 1), and, for s > 2, 

^z,i,o,i = ' uu ' f° r 3 < j < s + 1 (5.4) 

That is, for general s G I s , 

2 

Z(C n , q, s, v, w) = Y} X ^T + ( s - l )( vw T + (<? " s - !K ■ ( 5 - 5 ) 

It is readily checked that this expression for Z(C n , q, s, v, w) (i) reduces to the zero-field Potts 
model partition function 

Z(C n , q, v) = (q + v) n + (q- l)v n (5.6) 

for s — or w — 1, (ii) satisfies the general symmetry property f)2.13p . and (iii) reduces to 
w n Z(C n ,q,v) for s = q, in agreement with Eq. (12 . 181) . 

As was noted briefly at the beginning of the paper, our result f)5.5p shows a qualitative 
difference between the case s = 1 considered previously 0] and the more general set of cases 
with s > 2 in the interval I s , namely the fact that the third term, (s — l)(vw) n , is absent 
for s = 1 but is present for other values of s E I s . By the s <-> q — s symmetry inherent in 
Eqs. (12.41) and (I2.13p . this also means that another term vanishes identically for s = q — 1 
but is present for other values of s G I s ; this is the last term in (15.51) . (q — s — l)v n . It is 
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interesting to observe that the symmetry ( 12. 13ft applies not just to the total Z(C n , q, s, v, w), 
but also to parts of this function. Specifically, under the replacement s — > q — s, one sees 
that (i) the sum of the last two terms in ( 15. 51) . (s — l)(vw) n + (q — s — l)v n , transforms 



into w n [(q — s — l)v n + (s — ljfww" 1 )"] and (ii) the first two terms, ^2j = \[^z,i,a,j(q, s,v , w)] n 
transform into w n Ylj=i[^z,i,o,j(Q, Q ~ s,f,w _1 )] n , so that each of these parts, (i) and (ii), 
individually satisfies the symmetry (12.13p . 

We exhibit Z(C n ,q,s,v,w) for 2 < n < 4 below. To keep the equations as compact as 
possible, we write the coefficients of the terms of maximal degree in w and of degree in to 
in terms of zero-field partition functions using the general results ( 12.351) and (12.311) . We find 

Z(C 2 ,q, s,v,w) = Z(C 2 , s,v)w 2 + 2s(q - s)w + Z(C 2 ,q - s,v) 

= q 2 + [2t + v(v + 2)}q + t t + v(v + 2)(w + 1) (5.7) 

Z(C 3 ,q,s,v,w) = Z(C 3 , s,v)w 3 + 3s(q - s)(s + v)w 2 

+ 3s(q — s)(q — s + v)w + Z(C 3 , q — s,v) (5.8) 

Z(C4, q, s, v, w) = Z(Ci, s,v)w 4 + 4s (q — s)(s + v) 2 w 3 

+ 2s(q - s) 3(q 2 - s 2 - (q - s) 2 ) + 4qv + Av' 

+ Aq(q - s)(q - s + vfw + Z(C 4 , q - s,v) . (5.9) 

As usual, one obtains the Ph(C n ,q,s,w) for each n by setting v = — 1 in Z(C n ,q,s,v,w). 
For s — 1, the parts of Z{C n , q, 1, v, w) were given in Ref. 0] and Ph(C n , q, 1, w) was given 
inRef. Q- 

In the context of weighted-set coloring, so that s is an integer in the interval I s , it was 
noted [2j that for s — 1, Ph(C n , q, 1, w) contains the factor (q — 1). It was also noted that if 
n is odd, say n = 2m + 1 with m = 1,2,.. then Ph(C2 m +i, q,l,w) also contains the factor 
(q — 2), so that for odd n, P/i(C 2m +i, 2, 1, w) = 0]. The fact that if n = 2m + 1 is odd and 
q = 2, then one cannot perform a proper g-coloring of C n , so that -P/i(C 2m+ i, 2, s, w) = 0, is 
independent of both w and s. However, in contrast to the s = 1 case, where the vanishing of 
Ph{C2m+ii ?, 1, w) for 5 = 2 occurred via a factor of (q — 2), this is not the case for general 
s. Instead, Ph(C2 m +i, q, s, w) is such that if one evaluates it at q — 2, there is a factor which 
is a polynomial in s that implicitly but necessarily vanishes. This vanishing occurs because 



w 2 
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of the implicit restriction on the values of s, namely that s is an integer in the interval 
< s < q. For example, consider Ph(C^,q,s,w), given above as the v — — 1 special case of 
Eq. ( 15.81) . Evaluating this at q — 2, we obtain 

Ph(C 3 , 2, s, w) = s(s - l)(s - 2)(w - l) 3 . (5.10) 

This implicitly vanishes for any of the allowed (integer) values of s in the set I s because for 
q = 2, s can only take on the values 0, 1, or 2 in this set. The same type of mechanism 
is responsible for the vanishing of Ph(C2 m -\-i,q, s,w) at q = 2 for higher values of m. For 
s = 2, if n > 3 is odd, then our results show that Ph(C n , q, s, w) contains the factor (q — 2). 
If s > 3, then, in general, Ph(C n , q, s, w) does not have such overall factors. The reason for 
this is that s is bounded above by q, so that if s > 3, then q > 3. This is equal to x(C n ) 
for n odd and greater than x(C n ) for n even (cf. Eq. (15. ip ). so that Ph(C n ,q, s,w) is not 
forced to vanish the way it is for the cases of s — 0, 1, 2 when q can be less then x{C n )- 

We have mentioned above that our general structural formulas for Z(G,q,s,v,w) and 
Ph(G,q, s,w) with G a cyclic strip graph simplify considerably when s = q, as required 
by the general relations (12.181) and (ll.5p . It is interesting to see how this occurs for this 
family of circuit graphs, G = C n . There are Nzh,L y ,s = s + 2 terms \z,L y ,d,j(q, s, v , w) 
with L y — 1 whose n'th powers occur in Eq. (15.21) . Of these, nzh,i,o,s — s + 1 multiply 
the coefficient = 1 and the remaining term, [\z,i,i(q, s, v, w)] n = v n , multiplies the 
coefficient C« = q - s — 1. From the general relation (12.181) . we know that for s = q, 
Z(C n ,q,q,v,w) = w n Z(C n ,q,v), where Z(C n ,q,v) is the zero-field Potts model partition 
function for the circuit graph C n , given by Eq. (15. 6p . Hence, we can deduce that of the s + 1 
terms \z,i,o,j{Q, s > v , w ), (i) one becomes equal to w(q + v); (ii) a second becomes equal to v , 
and (iii) the remaining s — 1 = q — 1 terms become equal to wv. Using the fact that s = q, 
so that = —1, we then have the reduction 



s+1 

Z(C n ,q,q,v,w) = ^2[\z,i,o,j(q,q,v,w)] n + c 



w n 



w{q + v)] n + v n + (q- l){wv) n - v r 
(q + v) n + (q-l)v r 

= w n Z(C n ,q,v) . (5.11) 

Thus, the identity (12.181) (a special case of the symmetry (12.131) ) is realized via a "trans- 
migration" process in which one or more \z,L y ,dj(q, s, v , w)'s for a given d become equal or 
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proportional to (respectively, one or more) \z,L v ,d',j(q, s, v, w)'s for a different d! and hence 
their m'th powers can be regrouped with the latter, thereby changing the effective coeffi- 
cients that multiply the [Xz,L y ,d',j{q, s, v , w)] m, s , where here m = n. This transmigration 
process is a general one and occurs also for higher values of strip width L y . This process is 
also the mechanism whereby the identities (12. 13[) and (12.141) are satisfied. Thus, the results 
above for n Z h{L y , d, s), n Ph (L y , d, s), N Z h,L y ,s, and N PKLy , s apply for generic values of s, but 
there are simplifications, involving this type of transmigration, for special values of s. 

VI. SOME PROPERTIES OF THE ZEROS OF Ph(G,q,s,w) 

One can solve for the zeros of Z(G,q, s,v,w) as functions of any of the four variables 
q, s, v and w with the other three held fixed, and similarly, one can solve for the zeros of 
Ph(G, q, s, w) as a function of any of the three variables q, s, and w with the other two 
held fixed. Here it is understood that one uses Eq. (12.61) . with these four variables each 
generalized to lie in C. 

In general, some zeros in w and s can have unbounded magnitudes as a function of the 
other variables. The underlying reason for this can be seen at an algebraic level as the fact 
that the coefficient of the highest-degree term in this variable in Z(G, q, s, v, w) can vanish at 
some values of w and/or s. Related to this are the facts that (i) at the special values w = 1 
and w = 0, Z(G,q, s,v,w) loses its dependence on s, and (ii) at the special value s = 0, 
Z(G,q, s,v,w) loses its dependence on w. As a consequence, the zeros in these variables 
move off to infinity. Similar statements apply for Ph(G, q, s, w). We proceed to discuss some 
properties of these zeros. 

A. Zeros of Z(G,q, s,v,w) and Ph(G,q, s,w) in q 

Here we consider the zeros of Z(G, q, s, v, w) and Ph(G, q, w) in q, as a function of s and 
w, for some graphs G. Since the maximal degree of Z(G, q, s, v, w) and Ph(G, q, s, w) in the 
variable q is n(G), each of these polynomials has this number of zeros in the variable q. In 
general, since Z(G, q, s, v, 1) = Z(G, q, v) while Z(G, q, s, v, 0) = Z(G, q — s, v), it follows that 
as w decreases from 1 to 0, there is an overall shift to the right in the zeros of Z(G, q, s, v, w) 
in the q plane by s units, and this holds, in particular, for the case v — — 1 that yields 
Ph(G,q, s,w). This shift is illustrated by some simple examples. For G = L\ = K\ = N±, 
Z(Li, q, s, v , w) = at 

q Llz = -t = s(l - w) . (6.1) 
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This increases from to s as w decreases from 1 to in the DFSCP interval and decreases 
from through negative values as w increases above 1 in the FSCP interval. For L 2 = K 2 , 
Z(L 2 , q, s, v, w) vanishes at 

1 r 



qi2 



—v + 2s(l — w) ± y/v\v — Asw(w — 1)] 



(6.2) 



where j — 1, 2 for the ± sign. As w decreases from 1 to 0, the root qi2z,\ increases from to 
s, while the root qi2z,2 increases from — v to s—v. The zeros of graph-coloring polynomials in 
q satisfy certain boundedness properties for fixed nonzero magnetic field 2jJ. However, the 
zeros of Z(G, q, s, v, w) and Ph(G, q, s, w) in q are, in general, unbounded as — ^ oo. This 
is already evident in the simplest case of a single vertex, for which the zero of Z(L\, q, s, v, w) 
at q — s(l — w) has a magnitude that, for s ^ 0, goes to infinity as \w\ — > oo. 



B. Zeros of Z(G,q, s,v,w) and Ph(G,q, s,w) in s 

One may also study the zeros of Z(G,q, s,v,w) in each of the three variables s, v, and 
w with the other two (and q) held fixed, and the zeros of Ph(G, q, s, w) in each of the two 
variables s and w with the other (and q) held fixed. In general, these zeros are unbounded 
in magnitude even when the other variables vary over a finite range. The reason for the 
divergences in these zeros of Z(G,q, s,v,w) is the fact that the coefficient of the term in 
Z(G, q, s, v, w) of highest degree in the given variable can vanish as one changes other vari- 
ables, and similarly with Ph(G, q, s, w). Again, this can be illustrated with simple examples. 

We begin our discussion with the zeros in s, where it is understood that this variable is 
formally extended from the integers in the interval I s to the complex numbers. The Potts 
model partition function for the single-vertex graph L±, Z(Li, q, s, v , w), vanishes at 

1 

s = . 

1—W 

For q ^ 0, this diverges as w — > 1. For L 2 , Z(L 2 , q, s, v, w) has zeros in s at 



(6.3) 



-[2q + v(w + 1)] ± Jv[v(w + l) 2 + Aqw] , . 

SL ^ = 20^1) (6 ' 4) 

where j = 1,2 correspond to the ± signs, respectively. We find that 



jq + v) ± y/v{q + v) 

iz, 3 : asw^l (6.5) 

w — 1 

so that as w — > 1, the magnitudes |si 22J | are, in general, unbounded. These results for the 
full Potts model partition function apply, a fortiori, to the special case v = — 1 that defines 
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the weighted-set chromatic polynomial, Ph{L 2 , q, s, w). Above, we have explained the origin 
of this type of divergence as being due to the fact that the coefficient of the highest-power 
term of Z(G, q, s, v, w) and Ph(G, q, s, w) in the given variable, here s, vanishes. This occurs 
at w = 1. Closely related to this, at w = 1, all dependence on s in Z(G, q, s, v, w) disappears, 
so it is understandable that the zeros in s would disappear by moving off to infinity in this 
limit. 



C. Zeros of Z(G,q, s,v,w) and Ph(G,q, s,w) in w 

Here we comment on zeros of Z(G, q, s, v, w) and Ph{G, q, s, w) in w. First, we note that 
the symmetry (12. 13[) implies that if one replaces s by q — s, then the zeros of Z(G, q, s, v, w) 
and Ph(G, q, s, w) in w away from the origin map into their inverses. In particular, if q — 2s, 
then Eq. (I2.13P reads Z(G, 2s, s, v , w) = w n Z(G, 2s, s, v , u> _1 ), so that the zeros of Z away 
from the origin in the w plane form a set that is invariant under inversion. For G = L\, 
Z(Li, q, s, v, w) vanishes at 

wliz = 1 ■ (6.6) 

s 

From Eq. (13. 5J) , we find that Z(L 2 , q, s, v , w) = for 

s{s -q)± y/s(s - q)v(q + v) 

w L2z ,j = -7— — v (6.7) 

s(s + v) 

where j = 1,2 correspond to the ± signs, respectively. The right-hand sides of Eqs. (16.61) 
and (16.71) both diverge as s — > 0. In particular, 



^2,,i~±A/ V q)iq + V) ass^O (6.8) 
\ sv 

for j = 1,2. The fact that, in general, the magnitudes of the zero Wlu of Z{L\, q, s, v , w) 
and the zeros WL2z,j, j = 1,2, of Z(L 2 ,q, s,v,w) diverge as s — > is again understandable, 
since for a graph G, if s = 0, then Z{G, q, s, v, w) reduces to Z(G, q, v), with no dependence 
on w. Hence, in the limit s — > 0, it is natural that the zeros in w disappear by moving off to 
infinity. Finally, as s — > —v, we find that 

u) L 2z,2 -> 2v as s -> -v (6.9) 

while wl2z,i is unbounded: 

2(q + v) 

WL2z,i ~ as s — > —v . (6.10) 

s + v 
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D. Zeros of Z(G, q, s,v,w) in v 



To illustrate the calculation of the zeros of Z(G,q,s,v,w) in v, we again use our simple 
example graph, L 2 . We find the zero of Z(L 2 , q, s, v, w) in v occurs at 

[q + s(w-l)] 2 



q + s(w — l)(w + 1) 
This has an unbounded magnitude if q + s(w 2 — 1) = 0, i.e., 



(6.11) 



(6.12) 



1 — w 2 

This divergence in the magnitude of the right-hand side of (16.111) does not directly affect the 
weighted-set proper vertex coloring of this L 2 graph because in the DFSCP region < w < 1, 
the condition f 1 6 . X 2 f) implies that s > q, outside the actual coloring interval I s , and in the 
FSCP region w > 1, it implies that s is negative, again outside of this interval I s . 



VII. QUANTITIES DEFINED IN THE LIMIT n(G) -» oo 

A. / and <1> Functions 

Let us consider families of graphs G m that can be built up recursively, such as lattice 
strips. For such graphs, the m + 1 member of the family is obtained from the m th member 
by (possibly cutting and) gluing in a given subgraph. For example, for the graph C m , one 
cuts the circuit at some vertex and inserts another edge and vertex to get C m+ i, and so 
forth. Generalizing this, another example is a strip of the square lattice with transverse 
width L y , length L x = m, and periodic longitudinal boundary conditions, which we denote 
as sq(L y x L x ,cyc). In this case, the number of vertices is n = L y L x . Generically, the 
number of vertices of a recursive graph G m is of the form n(G m ) = am + b, where a and b are 
(integer) constants, so that the limit n — > oo is equivalent to the limit m — > oo. We denote 
the formal n — > oo limit of such graphs as {G} = linv^ooGm. In the present context, this 
n — > oo limit corresponds to the limit of infinite length for a strip graph of fixed width and 
some prescribed boundary conditions. Correspondingly, one can define the free energy per 
vertex as follows (with the subscript m suppressed in the notation) 

f({G}, s, v, w) = lim n~ l \n[Z(G, q, s, v, w)] (7.1) 

n— *oo 

and a function <3>({G}, q, s, w), 

$({G}, q, s, w) = lim [Ph{G, q, s, w)} 1/n . (7.2) 
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As before (cf. Eq. (1.9) of [22J and Eq. (2.8) of [23f| ) , one must take account of a noncommu- 
tativity of limits that can occur, namely the fact that for certain special values of q, denoted 
{q sp }, the limits n — > oo and q — > q sp do not commute: 

lim lim Z{G,q,s,v,w) 1/n ^ lim lim Z{G, q, s, v, w) 1/n (7.3) 

n-*oo q-^q sp q—*q sp n— >oo 

and the analogous formulas for the v = —1 case which defines &{{G},q, s,w). For further 
details, we refer the reader to our previous discussions of this j2[|22,|23|. An explicit example 
is provided by our result for Z(C n , q, s, v, w) in Eq. (15.51) ; if one sets q = s + 1 first before 
taking n — > oo, then the last term drops out, while if one takes n — > oo first with + 
then, since lim n ^ 00 [(g — s — l)v n ] 1 / n = v, the last term may remain in /. We see also an 
additional type of noncommutativity that is present, namely that if we extend s from an 
integer in I s to a real (or complex) variable, then for a set of special values of s, denoted 

lim lim Z{G,q,s,v,w) 1/n + lim lim Z{G, q, s, v, w) 1/n . (7.4) 

n — >oo s — >s sp s — >s sp n — >oo 

Again, our result for Z(C n , q, s, v, w) in Eq. (15.51) provides an illustration of this; if one sets 
s — 1 first before taking n — > oo, then the third term, (s — l)(vw) n , drops out, while if 
one takes n — ► oo first with s ^ 1, and then sets s = 1, it follows that, since lim n _ +00 [(s — 
l)(vw) n ] 1 / n = vw, the resultant limiting term vw may remain in /. Similarly, if one sets 
s = q — 1 before taking n — > oo, then the last term, (q — s — l)v n drops out, while if 
one takes n —>■ oo first with s ^ 9 — 1, and then sets s = q — 1, it follows that since 
lim n _ +00 [(g — s — l)v = v, the resultant limiting term v may remain in /. In our analysis 
of the n — ► oo limit for recursive families of graphs, unless otherwise indicated, we shall 
choose the order of limits in which we fix s first and then take n — > oo. 

The function Ph({G},q, s,w) generalizes the ground state degeneracy of the zero-field, 
zero-temperature Potts antiferromagnet, W({G},q) = lim^oo P(G, q) 1 ^. Thus (with 
care taken concerning the above-mentioned noncommutativity of limits), <&({G},q,s,l) = 
$({G}, q, 0, w) = W({G},q). In the case of the zero-field, zero-temperature Potts antifer- 
romagnet, the associated configurational entropy per vertex (which is thus the ground-state 
entropy per site) {G} is S — ks^TiW. The third law of thermodynamics states that the en- 
tropy per site S should vanish as the temperature goes to zero. However, there are a number 
of exceptions to this law. For instance, for the zero-field g-state Potts antiferromagnet on a 
square lattice, an elementary argument yields the lower bound S/Ub > (1/2) ln(g — 1), which 
is nonzero for q > 3. 

In the present case of the weighted-set chromatic polynomial, let us consider first the 
FSCP interval w > 1 and assume that q > 2, so that a proper g-coloring can be performed 
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for a bipartite graph. We also assume that s is an integer for this discussion and lies in the 
interval 2 < s < q. Applying our lower bound (I2.78P to a particular bipartite graph, namely 
the strip graph of the square (sq) lattice with width L y vertices and even length L x vertices, 
denoted sq(L y x L x ), we have (with n\ = n<i = L x L y /2) the lower bound 

Ph(sq(L y x L x ), q, s, w) > s(s - l) n/2 w n . (7.5) 

Taking L x — > oo with L y fixed, we thus obtain the lower bound Q(sq(L y x oo), q, s, w) > 
wy/s — 1. Hence, in this limit, the entropy is bounded below by 

S(sq(L y x oo), q,s,w) > \nw + ^ ln(s - 1) (FSCP case) . (7.6) 

For s > 3 (which implies q > 3 also), this ground state entropy is nonzero. In the DFSCP 
interval < w < 1, with 1 < s < q — 3, applying our lower bound ( 12.791) to the same limit 
of this strip graph, we have &(sq(L y x oo), q, s, w) > y/q — s — 1, so that 

S{sq{L y x oo), q,s,w) > -ln(q-s- 1) (DFSCP case) . (7.7) 
For the given range, q > s + 3, this entropy is again nonzero. 



B. Example of Calculation of / and for a Family of Graphs 

We illustrate the calculation of the functions / and $ for the n — > oo limit of the circuit 
graph C n , or equivalently, the one-dimensional lattice with periodic boundary conditions. We 
shall take s to be a fixed integer in the interval < s < q for this analysis. For a given set of 
values of q, s, v, and w, the functional form of / is determined by the term \z,i,o,j{q, s, v, w) 
in Eq. (15.51) with the largest magnitude. For fixed s, v, and w and sufficiently large real q, 
this is Az.1,0,1- Following our nomenclature in earlier work for w — 1, we denote this region 
as region R\. As for the zero- field case, / and $ in this region are the same for the n — ► oo 
limit of the line graph L n and the circuit graph C n . We denote these limits as {L} and {C}. 
We thus have 

f({L}, q, s, v, w) = f({C}, q, s, v, w) = m[A^i, 0) i(g, s, v, w)) (7.8) 

and 

$({L}, q, s, w) = $({C}, q, s, w) = ln[A PM , ,i(g, s, w)} (7.9) 

where \z,i,o,j(q, s,v,w), j = 1,2, were given for this family of graphs in Eq. (15. 3p and 

ApM,o,j(g, s > w ) = x z,i,o,j(q, s, -1, w). 
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It is of interest to consider how $>(G,q,s,w) for {G} = {L} (or equivalently {G} = {C} 
in region Ri) behaves for certain special cases or limits of its variables. For example, in the 
limit where the weighting is removed, i.e., for w — > 1, <3>({L},g, s,w) has the Taylor series 
expansion 

^/rn \ -, s(q -l)(w-l) s(q - l)(q - s)(w - l) 2 , s \ 

$({L},g,s,w) = g-l + — - - „ -+0 [{w-1 3 asw -> 1 . 

g g^ V / 

(7.10) 

For s = 1, it was shown in Ref. j^J that the asymptotic behavior of $({L}, g, s, w) for large 
\w\ is 



$({L},g,s,w) ~ v/(g- l)w 
In contrast, for s ^ 1, we find 



1 



for s = 1 and Iwl — > oo . 



(7.11) 



$({L},g,s,w) ~ (s - l)w + - — S - — - + 



s(q — s) + q — 1 
2(s-l) 



1(7 



(7.12) 



As is evident, the large-|w| behavior is different, depending on whether or not s = 1. For 
|g| — > oo, we obtain the asymptotic expansion 



$({L}, q, s, w) ~ q + s(w - 1) - 1 - 



sw(w — 1) 
Q 



Oil 



(7.13) 



Extending s from an integer in the interval I s to a real (or complex) number, it is of interest 
to determine the limiting behavior of $({L}, q, s, w) as s — > and s — > g. We calculate 



$({L},g,s,w) = g-l + 



(w - 1)(g - 1)g +0( g 2 ) as, 
ty + g - 1 







(7.14) 



and 



<S>({L},q,s,w) = w(q-l)- W{W } )iq ^ S) +0((q-s) 2 ) as S - g (7.15) 

w(g — 1) + 1 

From our study of $({G}, g, s, w) for the n — > oo limits of various families of graphs, we 
have observed several generic properties. Let us consider families of strip graphs G of regular 
lattices A. Let the chromatic number of the lattice A be denoted as x(A) and assume that 
g > x(A). A technical assumption is that the n — > oo limit of the lattice strip graphs is 
taken in a manner such that for each G, x{G) = x(A). (For example, for square-lattice strip 
graphs with periodic longitudinal boundary, this means taking the length to be even.) Then 
for fixed s G I s we have observed that (i) for fixed w > 0, $({G}, g, s, w) is a monotonically 



39 



increasing function of q and (ii) for fixed q, $({G}, q, s, w) is a monotonically increasing 
function of w for w > 0. One can also analyze the behavior of $({G}, q, s, w) as a function 
of s for fixed q and w, but because of the noncommutativity (17.41) . one must take care in 
specifying the order of limits used in defining this function. The order that we take is first 
to set s to a given value in I s and then to take n —>■ oo. As is evident in the definition of 
I s , it is understood that s < q. We find that if w > 1, then $({G}, q, s, w) is an increasing 
function of s G I s , while if < w < 1, then q, s, w) is a decreasing function of s G I s . 

Focusing on the circuit graph, one sees that for values of the variables such that another 
term A in Eq. (15. 5p becomes dominant, there is a non- analytic change in / and $. As we have 
discussed earlier, this is also associated with a locus, denoted generically B, that comprises 
the accumulation set of zeros of the respective function, / and For definiteness, we analyze 
the locus B in the q plane (denoted B q ) for $({C}, q, s, w ), with s and w fixed. It may be 
recalled that for the n — > oo limit of the unweighted chromatic polynomial, B q is the unit 



circle \q — 1| = 1, so q c = 2 in that case 22|, |26j|. For the weighted-set chromatic polynomial 



Ph(C n , q, s, w), the accumulation locus B depends on the value of s and on whether w is in 
the DFSCP interval < w < 1 or the FSCP interval w > 1. We consider here the DFSCP 
interval, since as w decreases from 1 to 0, Ph(C n , q, s, w) interpolates between P(G, q) and 
P(G,q — s). We also generally take s ^ 0, since for s = Ph(C n ,q, s,w) reduces to the 
well-studied unweighted chromatic polynomial P(C n , q). (However, our results subsume this 
s = case.) 

In this DFSCP interval, for s = 1 or s = 2, the lower boundary of the region R\ on the 
real axis, denoted q c , is determined by the equality in magnitude 



I^PMAil — |<\p/i,i,i| — 1 ' (7-16) 



which yields the result 



s ( 1 — w ) 

q c = 2 + -4 '- for {G} = {C} and s = 1 or s = 2 and < w < 1 . (7.17) 

1 + w 

For < w < 1, this value of g c is greater than the value g c = 2 for the unweighted chromatic 
polynomial. Furthermore, as is evident from Eq. ( 17.171) . q c is a monotonically increasing 
function of s for fixed u> in this DFSCP interval and a monotonically decreasing function of 
w for the fixed values of s given above. As w decreases from 1 to 0, q c increases continuously 
from 2 to 2 + s. In contrast, the left-hand part of the boundary changes discontinuously; as 
w decreases by an arbitrarily small amount below 1, the point on the left where B q crosses 
the real q axis jumps discontinuously from q = to q = s. These results are in accord with 
the fact that for w = 0, B q is the locus of solutions to the equation \q — (1 + s)\ — 1, i.e., 
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the unit circle in the q plane centered at the point q = 1 + s, crossing the real axis on the 
left at q = s and on the right at q = 2 + s. The locus B g separates the q plane into two 
regions. We label the regions outside and inside the closed curve B as R\ (noted before) and 
i?2, respectively. 

For s > 2, there is a change in the locus B q , because the condition of degeneracy of leading 
A's is different; rather than Eq. (I7.16p . it takes the form 

|<\pMAil = I^PMA2| • (7-18) 

This entails the condition that (i) 

q- s - 1 + w(s - 1) = , (7.19) 

so that Xph,i,o,i = —Xph,i,o,2, and the condition that (ii) lAp^i^il = (Ap^o^l > 1, so that 
these A's are dominant. Substituting for q from Eq. ( 17. 19ft . we find that condition (ii) is 
satisfied in the relevant range of w for 

— - < w < 1 . (7.20) 
s — 1 

This interval is nonvanishing if s > 2. Thus, for s > 2, provided that conditions (i) and (ii) 
are satisfied, B q has the form of an open self-conjugate arc that crosses the real axis at the 
point given by Eq. (I7.19p . so that in this case, 

q c = 8 + 1 - w (s - 1) for {G} = {C} and s > 2 and < w < 1 . (7.21) 

s — 1 

This self-conjugate arc is concave to the left and ends at the arc endpoints given by where 
the (v — — 1 evaluation of the) expression in the square root of Eq. (I5.3P vanishes, namely 

q e , q* = (s + - w) ± 2iy/sw(l - w) . (7.22) 

This locus does not separate the q plane into different regions. It is straightforward to carry 
out a similar analysis of B q for the FSCP regime w > 1. 

VIII. RELATED TOPICS 

As our results show, one finds a number of intriguing features in the study of weighted-set 
vertex coloring of graphs. There are many further directions of research in this general area. 
One could, for example, use the methods presented here to calculate Z(G,q, s,v,w) and 
Ph(G, q, s, w) for other individual graphs and families of graphs. One could also investigate 
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further the zeros of these functions in various variables and their accumulation sets B for 
recursive graphs in the limit of infinitely many vertices. It would, moreover, be worthwhile 
to study connections with weighted loop models [271 ] . One could also investigate a different 
but related type of graph coloring problem in which the set of colors that one chooses from to 
assign to each vertex depends on the vertex. The unweighted case is called the list coloring 
problem in graph theory [28(, and it would be useful to study the weighted-set generalization 
of list coloring. We are pursuing these studies. 



IX. CONCLUSIONS 



In this paper we have studied the weighted-set graph coloring problems, in which one 
assigns q colors to the vertices of a graph such that adjacent vertices have different colors, 
with a vertex weighting w that either disfavors or favors a given set of s colors. In particular, 
we have analyzed an associated weighted-set chromatic polynomial Ph(G, q, s, w) and have 
also related this to a corresponding Potts model partition function with external magnetic 
fields, Z(G,q,s,v,w). These functions exhibit a wealth of interesting properties. We have 
proved various general results on these and illustrated them for particular graphs and families 
of graphs. 
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